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A  mixed  formulation  for  calculating  static  equilibrium  and  stability 
eigenvalues  of  nonuniform  rotor  blades  in  hover  is  presented.  The  static 
equilibrium  equations  are  nonlinear  and  are  solved  by  an  accurate  and  effi¬ 
cient  collocation  method.  The  linearized  perturbation  equations  are  solved 
by  a  one-step,  second-order  integration  scheme.  The  numerical  results 
correlate  very  well  with  published  results  from  a  nearly  identical  stabil¬ 
ity  analysis  based  on  a  displacement  formulation.  Slight  differences  in 
the  results  are  traced  to  terms  in  the  equations  that  relate  moments  to 
derivatives  of  rotations.  With  the  present  ordering  scheme,  in  which  terms 
of  the  order  of  squares  of  rotations  are  neglected  with  respect  to  unity, 
it  is  not  possible  to  achieve  completely  equivalent  models  based  on  mixed 
and  displacement  formulations.  A  study  of  the  one-step  methods  reveals  that 
a  second-order  Taylor  expansion  is  necessary  to  achieve  good  convergence 
for  nonuniform  rotating  blades.  Numerical  results  for  a  hypothetical  non- 
uniform  blade,  including  the  nonlinear  static  equilibrium  solution,  were 
obtained  with  no  more  effort  or  computer  time  than  that  required  for  a 
uniform  blade  with  the  present  analysis. 


1 .  Introduction 


It  has  been  found  that  nonlinearities  in  rotor-blade  equations  affect 
blade  stability  [1-4]  —  especially  stability  of  coupled  flap,  lead-lag,  and 
torsion  degrees  of  freedom  [1].  In  these  references  it  was  found  that 
essential  nonlinear  effects  could  be  retained  by  perturbing  the  nonlinear 
equations  of  motion  about  the  static  equilibrium  condition  and  solving  for 
the  eigenvalues  of  the  linearized  perturbation  equations.  Coefficients  of 
the  linearized  perturbation  equations  then  depend  on  the  solution  of  the 
nonlinear  static  equilibrium  equations. 


This  paper  presents  a  method  for  stability  analysis  of  nonuniform 
rotating  blades  with  aerodynamic  loading  that  utilizes  the  solution  for  the 
nonlinear  static  equilibrium  equations  in  the  linearized  eigenvalue  problem. 
Methods  described  in  the  literature  have  been  limited  to  solution  of 
various  restricted  versions  of  this  problem  [4-9].  These  include  a  modal 
approach  [4],  an  integrating  matrix  method  [5],  a  Myklestad  method  [6], 
a  Ritz  finite-element  method  [7],  a  modal  approach  based  on  a  mixed  varia¬ 
tional  principle  [8],  a  transmission  matrix  method  [9],  and  a  Galerkin 
finite-element  method  [10]. 
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In  section  2  of  this  paper,  as  in  [1],  the  governing  equations  of  motion 
are  written  as  two  sets  of  equations.  One  is  a  set  of  nonlinear  ordinary 
differential  equations  that  governs  the  static  equilibrium  condition.  The 
other  is  a  set  of  linear  ordinary  differential  equations  with  an  unknown 
eigenvalue  that  governs  the  dynamic  behavior  of  small  perturbation  motions 
about  equilibrium.  These  differential  equations  are  written  in  a  mixed, 
spatial-derivative  form,  unlike  the  displacement  formulation  of  [1]. 

Several  differences  between  mixed  and  displacement  formulations  are  then 
discussed  in  the  context  of  the  present  analytical  task. 

Then,  in  section  3,  solution  of  the  equations  is  discussed.  A  differ¬ 
ent  technique  is  used  to  solve  for  the  nonlinear  static  equilibrium  from 
that  used  to  solve  for  the  linear  stability  eigenvalues.  The  nonlinear 
static  equilibrium  is  solved  by  a  collocation  method  for  a  mixed-order  sys¬ 
tem  of  boundary  value  equations  [11].  The  software  for  this  method  is 
available  in  a  program  called  COLSYS  [12].  For  the  linearized  perturbation 
equations,  several  methods  of  solving  for  the  eigenvalues  are  discussed. 

These  include  COLSYS  [12],  a  generalized  eigenvalue  approach  [13],  a  Block- 
Stodola  technique  [14],  and  one-step  numerical  integration  techniques  of 
various  orders.  For  the  present  formulation  the  one-step  methods,  which 
appear  to  be  the  most  promising,  are  used. 

Finally,  in  section  4,  numerical  results  are  presented  for  comparison 
with  published  data  and  for  study  of  convergence  properties  associated  with 
the  one-step  methods.  Numerical  results  are  also  presented  for  stability 
eigenvalues  of  a  hypothetical  nonuniform  blade. 

2.  The  Governing  Equations 

The  governing  equations  used  herein  were  essentially  derived  in  [15], 
in  which  the  Houbolt  and  Brooks  [16]  equations  are  extended  into  the  range 
of  geometric  nonlinear  behavior.  The  integro-partial  differential  equa¬ 
tions  are  transformed  into  ordinary  differential  equations  by  first  express¬ 
ing  the  vector  of  unknowns,  z,  as 

z(x,t)  -  z(x)  +  z(x,t)  (1) 

Here  z  is  the  static  equilibrium  part  and  z  is  a  small,  perturbed  part. 
The  equations  are  linearized  in  z  and  converted  to  an  eigenvalue  problem 
with  the  transformation 

z(x,t)  -  z(x)eXt  (2) 

where 

A  •  o  +  iw  (3) 

The  static  equilibriua  equations  (with  the  sign  conventions  and  nomenclature 
defined  in  [15]  are  given  in  appendix  A.  From  Equations  (Al)  through  (A10), 
and  (A17),  the  state  vector,  z,  for  the  nonlinear,  static  equilibrium 
behavior  is  written  as 

z  -  [Vy  Mz  v  l  V2  My  w  6  Mx  J  ij1  (4) 

and  Equations  (A13) ,  (A14) ,  and  (A15)  define  the  aerodynamic  loading  for 
the  hovering  flight  condition. 
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The  linearized  equations  governing  small  perturbation  motions  from  the 
static  equilibrium  state  defined  by  the  solution  for  the  equations  in  appen¬ 
dix  A  are  given  in  appendix  B.  The  state  variables,  z,  for  the  perturbed 
state  are  defined  in  Equations  (Bl)  through  (BIO),  and 
The  vector  z  is 

z  =  I V  M  v  l  VMw  $  M  4>  v*  w* 

L  y  z  z  y  x 

and  the  equations  may  be  represented  in  matrix  form  as 

=  (A  +  AB)z  (6) 

In  the  derivation  of  the  equations  in  appendices  A  and  B,  terms  of  0(e2) 
have  been  dropped  with  respect  to  unity  where  e  is  the  order  of  bending 
and  torsion  rotations,  £,  3,  and  <p.  The  ordering  scheme  outlined  in  [15] 
has  been  followed  as  closely  as  possible.  As  in  [15],  an  exception  has 
been  made  to  include  0(e2)  terms  in  the  torsion  equations  that  are  uncou¬ 
pled  from  the  bending  equations  but  are  known  to  influence  the  uncoupled 
torsion  frequency.  It  is  not  possible,  however,  to  be  completely  consistent 
in  ordering  schemes  in  either  a  mixed  formulation  for  the  differential 
equations  or  in  a  displacement  formulation.  Furthermore,  as  will  be  shown 
in  the  next  section,  it  is  not  always  possible  to  reach  complete  agreement 
between  the  mixed  and  displacement  formulations  in  the  decision  of  which 
terms  to  retain,  even  when  trying  to  follow  the  same  guidelines  on  neglect¬ 
ing  higher-order  terms. 

An  exception  to  the  guidelines  in  [15]  on  the  ordering  of  terms  has 
been  to  retain  complete  trigonometric  expressions  involving  the  variable  <f>. 
In  [15],  cos(0  +  <J>)  and  sin(6  +  <Ji)  were  approximated  by  the  appropriate 
expansions  in  terms  of  cos  0,  sin  0,  and  4> ,  neglecting  terms  of  second  and 
higher  powers  of  <j>.  In  this  paper  it  has  been  convenient  to  retain  the 
complete  trigonometric  expression  in  all  equations  except  the  equations 
relating  rotation  derivatives  and  moments  which  are  presented  below.  No 
significant  differences  in  the  results  of  [1]  and  this  paper  will  occur  as 
a  consequence  of  using  this  convenient  trigonometric  form. 

When  the  equations  are  formulated  as  a  system  of  first-order  differen¬ 
tial  equations  (i.e.,  a  mixed  formulation)  several  attractive  features 
become  apparent.  These  include  the  simplicity  of  the  form  of  the  governing 
differential  equations,  the  absence  of  the  derivatives  of  the  elastic  char¬ 
acteristics,  the  simplicity  of  applying  numerical  integration  techniques, 
and  of  handling  boundary  and  interface  conditions.  The  mixed  formulation 
was  particularly  convenient  in  this  application  since  the  integro- 
differential  terms  could  easily  be  rewritten  as  differential  equations 
[note  Eq.  (A17) ,  (B14) ,  and  (B19) ] . 

Concerning  retention  of  specific  terms,  the  main  difference  between  a 
mixed-order  formulation  and  displacement  formulations  involves  the  equations 
relating  derivatives  of  rotations  to  moments  which  are 


(B14)  through  (BI9). 

*  *  *  It 

<t>*  vx  uj  (5) 


aiC  -  a26  -ai  +  a^J 
(s')  a3;  -  ai6  -a 3  -  2ai$ 


a2 

ai 


”x 

M„ 


+ 


<eAVx/Z0> 


(7) 
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These  equations  lead  to  Equations  (A4)  and  (A8)  in  appendix  A  and  Equa¬ 
tions  (B4)  and  (B8)  in  appendix  B.  Reference  [1],  which  is  a  displacement 
formulation  neglects  the  M  contribution  (i.e.,  the  single-underlined  terms 
in  Equation  7)  in  a  consistent  manner  on  the  basis  that  torsional  moments  are 
at  least  one  order  of  magnitude  less  than  the  bending  moments.  References  [1] 
and  [15]  use  a  quasi-coordinate  as  the  torsion  variable  as  discussed  in  [17]. 
For  that  formulation,  the  torsion  moments  when  used  in  the  bending  equations 
are  written  as  integrals  of  applied  and  inertial  loads  and  are  thus  of 
0(e2).  In  the  bending  equations  they  are  multiplied  by  quantities  of  0(e) 
and  are  thus  negligible.  References  [18]  and  [19],  which  are  also  displace¬ 
ment  formulations,  use  Lagrangian  torsion  variables  and  in  a  consistent 
manner  retain  the  single-underlined  terms.  There,  the  torsion  moments  are 
written  in  terms  of  the  first  derivatives  of  4>  and  are  0(e)  quantities. 

When  they  are  multiplied  by  terms  of  0(e)  in  the  bending  equations,  they  are 
not  formally  negligible  and  thus  are  retained.  The  single-underlined  terms 
are  also  retained  in  the  mixed  formulation  presented  here,  although  they  could 
be  neglected  on  the  basis  of  the  arguments  in  [1] .  Results  are  presented  in 
a  later  section  which  document  the  effect  of  the  single-underlined  terms. 

Another  more  important  difference  between  mixed  and  displacement  for¬ 
mulations  is  illustrated  by  the  double-underlined  terms  of  Equation  (7). 

In  [1],  [15],  [18],  and  [19]  the  double-underlined  terms  were  retained  in 
deriving  the  fourth-order  governing  equations  for  the  lead-lag  and  flap¬ 
bending  behavior,  but  they  were  neglected  in  deriving  the  second-order 
torsion  equation.  These  terms  are  consistently  retained  or  neglected  in 
the  above  references  on  the  basis  of  essentially  the  same  ordering  scheme, 
and  the  retained  terms  result  in  a  symmetric  elastic  stiffness  matrix.  For 
the  present  mixed  formulation,  however,  there  is  no  reason,  a  priori,  to 
neglect  these  terms.  We  may  expect  these  terms  to  have  some  observable 
effect,  especially  when  there  is  significant  static  equilibrium  deformation, 
such  as  when  0  is  large.  This  effect  will  be  illustrated  in  the  next 
section. 

Previously,  we  discussed  some  advantages  of  the  mixed  formulation.  There 
are  also  some  drawbacks  to  use  of  mixed  formulations.  The  matrix  operator  for 
the  structural  terms  in  the  equations  no  longer  appears  in  symmetric  form,  and 
the  number  of  variables  in  the  state  vector  is  several  times  that  of  displace¬ 
ment  formulations.  For  some  applications,  however,  the  convergence  is  more 
rapid  for  mixed  formulations  than  for  displacement  formulations  [20  and  21], 
thereby  offsetting  somewhat  the  latter  disadvantage.  Further,  force  resultants 
are  obtained  with  the  same  degree  of  accuracy  as  the  displacements  in  the 
mixed  formulation.  The  disadvantage  of  the  nonsymraetric  form  of  the  structural 
operater  can  be  somewhat  offset  by  selection  of  solution  techniques  that  do 
not  make  use  of  matrix  symmetry  or  positive  definiteness.  Basically,  the 
choices  considered  here  were  whether  to  use  collocation  methods,  generalized 
eigensystem  approaches,  or  one-step  integration  methods. 

3.  Solution  of  the  Systems  of  Equations 

For  obtaining  the  nonlinear,  static  equilibrium  solution,  the  colloca¬ 
tion  method  [11]  was  used  since  the  software  [12]  was  readily  available 
and  since  the  method  ensures  a  high  accuracy.  The  program  is  capable  of 
treating  general  systems  of  nonlinear  multipoint  boundary  value  problems  up 
to  order  four  with  a  variety  of  options  available  to  the  user.  COLSYS  pro¬ 
duces  an  approximate  solution  to  a  user-specified  accuracy  using  a  polynomial 
spline  which  can  be  evaluated  at  any  point  within  the  domain  of  the  problem. 

The  Inertial,  geometric,  and  elastic  properties  of  the  blade  can  be  expressed 
as  functions  of  the  axial  coordinate  in  this  approach,  and  the  high-order. 
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spline-fit  method  leads  to  almost  as  high  a  degree  of  accuracy  at  noncollo¬ 
cation  points  as  at  the  collocation  points.  Further,  trigonometric  expres¬ 
sions  involving  state  variables,  such  as  sin(6  +  <j>),  can  be  expressed 
exactly  without  the  small  angle  assumptions  on  $  that  are  often  made. 

The  collocation  software  package,  COLSYS,  was  written  to  solve  a  set 
of  nonlinear  multipoint  boundary  value  problems.  As  such,  it  can  also  be 
applied  to  eigenvalue  problems  using  the  approach  outlined  in  [22].  After 
investigation  of  this  approach  for  the  eigenvalue  problem,  it  was  discarded 
for  various  reasons.  In  particular,  since  the  complex  roots  doubled  the 
number  of  required  governing  equations,  the  existing  code  could  no  longer 
accommodate  these  cases  without  modification.  A  further  limitation  is  that 
a  priori  estimates  of  the  desired  complex  eigenvectors  and  eigenvalues 
must  be  provided. 

Next,  the  generalized  eigenvalue  approach  was  investigated.  In  this 
approach.  Equation  (6)  is  discretized  using  a  finite-difference  approxima¬ 
tion.  The  one  chosen  here  involved  the  central  difference  approximation 
given  by: 


'i+i 


-  z , 


(h/2)<*  +  XB)1+l/!(,1  +  z1+i) 


(8) 


where  the  subscripts  denote  the  locations  at  which  the  variables  are 
evaluated. 

This  results  in  a  generalized  eigenvalue  problem  of  the  form 

•  (A  +  AB)Z  -  0  .  (9) 

The  dimensions  of  the  matrices  A  and  B  are  16k  *  16k,  where  16  is  the 
dimension  of  the  matrices  A  and  B  In  Equation  (6) ,  and  k  is  the  num¬ 
ber  of  segments  in  the  discretization.  In  preliminary  investigations  the 
subroutine  EIGZF  [13]  was  used  to  obtain  all  the  eigenvalues  for  this 
generalized  eigenvalue  problem.  It  was  found  that  as  many  as  100  segments 
would  be  required  to  obtain  3-place  accuracy  for  the  second  eigenvalue. 
Unfortunately,  the  lack  of  symmetry  prevents  the  banded  structure  of  the 
matrices  from  being  exploited;  thus,  a  storage  problem  is  created. 

One  alternate  approach  is  to  find  selected  eigenvalues.  The  Block- 
Stodola  method  [14]  uses  a  block  inverse-iteration  algorithm  to  find  the 
first  few  eigenvalues.  In  those  cases  which  were  tested  for  this  method, 
no  more  than  the  first  eigenvalue  could  be  computed  accurately.  Here, 
because  of  the  lack  of  symmetry,  there  is  no  underlying  variational  princi¬ 
ple  that  can  be  employed  as  in  [14].  Although  it  is  possible  to  symmetrize 
one  of  the  two  matrices  and  obtain  the  hypotheses  in  [14],  this  results  in 
a  significant  degree  of  fill-in.  As  a  consequence,  this  approach  was  not 
pursued. 

The  one-step  methods  using  transfer  matrices  similar  to  the  Myklestad 
method  seem  to  be  the  most  promising.  The  term  one-step  method  refers  to 
those  methods  which  depend  only  upon  station  i  to  obtain  solutions  at 
station  i  +  1.  Higher-order  terms  of  a  Taylor  expansion  can  be  used  with 
the  one-step  method  to  speed  convergence.  The  method  can  be  described  as 
follows.  A  Taylor  expansion  of  z(x^+j)  about  x^  is 

z<x  )  -  z(x^)  +  z'U^h  +  z" (x^)h2 /2  +  z'”(xi)h3/6  +  .  .  .  (10) 
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where 


h 


i+i 


-  x. 


Now,  in  Equation  (6) ,  let 


C  -  A  4-  XB 


(11) 


(12) 


A  first-order  method  can  be  obtained  from  the  first  two  terms  of  Equa¬ 
tion  (10),  yielding 


zi+x 


zi  + 


hCizi 


(13) 


where  Equation  (6)  is  used  to  eliminate  z'(xi).  Equation  (13)  has  0(h2) 
local  truncation  error.  A  second-order  method  can  be  similarly  obtained 
from  the  first  three  terms  in  Equation  (10)  yielding 


z 


141 


Zi  +  M  iZi  +  <h2/2Kcizi  +  ci2!> 


(14) 


which  has  0(h3)  local  truncation  error.  Here  Equation  (6)  is  used  to 
eliminate  z'(xj)  and  z"(x^).  Both  Equations  (13)  and  (14)  may  be  rewritten 
as  z1+1  ■  Tjzi,  where  is  a  suitably  chosen  matrix  evaluated  at  xi, 

so  that 


zi4-i  “  (I  +  hCi)zi  (15) 

z1+i  -  [1  +  hC1  +  (h2/2)(C^  4-  Cj)]2i  (16) 

The  form  of  T^  depends  on  the  order  of  the  method.  Hence,  zn  *  Tzo  where 
T  *  Tn-iTn_2 . . ,T0.  The  homogenous  boundary  conditions  at  Xq  and  xn  are  then 
used  to  reduce  matrix  T  to  a  smaller  (6  x  6)  matrix  T  whose  determinant 
det  (T)  is  a  polynomial  in  X.  The  desired  eigenvalues  correspond  to  those 
X *  8  that  produce  zero  determinant  for  f.  In  our  Implementation  of  the 
one-step  methods,  subroutine  ZANLYT  [13]  is  used  to  find  one  or  more  of  the 
complex  roots  of  the  real  golynomlal  det(T).  Subroutine  LEQTIC  [13]  is 
used  to  decompose  matrix  T  into  L  and  U  factors.  The  product  of  the 
pivots  produced  from  this  factorization  gives  the  determinant  of  T.  Finally, 
whenever  is  used,  it  is  approximated  by  the  backward  difference 

CI  *  <Ci  ‘  Ci-i)/h  (17) 

and  Co  is  assumed  to  be  zero.  In  any  one-step  method,  care  must  be  taken 
to  define  the  C'  terms  properly  at  those  points  where  cusps  or  discon¬ 
tinuities  in  the  axially  varying  properties  occur. 

If  there  are  axial  variations  in  inertial,  geometric,  and  elastic 
properties  or  in  the  tension,  the  term  C^z^  in  Equation  (16)  can  have  sig¬ 
nificant  Impact  on  the  speed  of  convergence.  Rotor  blades  have  all  of 
these  axial  variations.  Therefore,  the  method  illustrated  by  Equation  (16) 
was  used  in  this  paper.  As  shown  in  the  next  section  the  number  of  seg¬ 
ments  required  to  produce  accurate  results  increases  substantially  if  the 
Cjz^  term  is  omitted  from  Equation  (16). 
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4.  Numerical  Results  and  Discussion 


Numerical  results  were  obtained  using  an  IBM  360/67  computer  and  are 
presented  in  this  section  for  several  rotor-blade  configurations.  These 
results  are  intended  to  serve  primarily  three  purposes:  (1)  to  compare 
with  published  results  and  thus,  at  least  partially,  validate  the  present 
computer  program;  (2)  to  study  the  convergence  properties  of  one-step 
methods;  and  (3)  to  present  some  new  results  for  a  nonuniform  blade  which 
may  serve  as  a  reference  problem  for  future  analytical  studies.  The 
numerical  values  of  the  various  inertial,  elastic,  and  geometric  properties 
are  presented  in  Table  1  for  configurations  to  be  compared  with  results 
from  [1]  and  [9]  and  for  the  hypothetical  nonuniform  blade. 

4.1.  Comparison  with  published  results 

In  this  section  numerical  results  are  presented  to  compare  with  some 
published  data.  No  attempt  will  be  made  to  compare  with  all  of  the  many 
available  numerical  results  in  the  literature.  Instead,  we  will  focus  on 
results  from  an  in  vacuo  configuration  in  [9]  and  the  aeroelastic  stability 
results  of  [1],  Because  of  the  very  good  agreement  between  the  present 
results  and  those  published,  it  was  felt  that  a  presentation  of  the  results 
in  tabular  form  would  facilitate  comparison  with  published  data  and  under¬ 
standing  of  the  effects  of  the  underlined  terms  in  Equation  (7)  on  the  aero¬ 
elastic  stability  results  of  [1].  For  serving  the  latter  purpose  we 
designate  two  constants  kl  and  k2.  For  kl  (or  k2)  =  1,  the  single- (or 
double-)underlined  terms  in  Equation  (7)  are  included.  For  kl  (or  k2)  =  0, 
the  single-(or  double-)underlined  terms  in  Equation  (7)  are  deleted. 

As  pointed  out  above,  consistent  application  of  the  ordering  scheme  in 
a  displacement  formulation,  such  as  in  [1],  may  lead  to  discarding  terms 
that  may  not  necessarily  be  negligible  in  an  analogous  ordering  scheme  for 
a  mixed  formulation  as  in  this  paper.  In  order  to  match  results  of  [1],  it 
was  necessary  to  experiment  with  deleting  and  retaining  the  underlined  terms 
in  Equation  (7).  One  example  of  the  dilemmas  faced  in  trying  to  be  con¬ 
sistent  is  apparent  when  the  torsion  moment  Mj^  is  written  as  an  integral 
of  applied  and  inertial  loads.  In  this  case  it  is  clearly  0(e2),  as  noted 
in  [15]  and  explained  in  [17].  However,  when  Mx  is  written  in  terms  of 
the  derivatives  of  rotations,  it  contains  one  term,  which  is  0(c).  If  Mx 
is  regarded  as  0(e),  the  single-underlined  terms  must  remain  (kl  =  1);  if 
Mx  is  regarded  as  0(e2),  these  terms  should  be  neglected  for  the  sake  of 
consistency  (kl  =0).  A  second  example  involves  the  double-underlined 
terms.  When  the  moment  components  My  and  Mz  appear  in  the  bending  equa¬ 
tions  in  a  displacement  formulation,  they  appear  multiplied  by  0(1) 
quantities.  On  the  other  hand,  in  the  torsion  equation  they  are  multiplied 
by  0(e)  quantities.  Because  terms  up  through  0(e2)  are  retained  in  all 
equations  in  [1],  it  is  consistent  to  take  only  the  dominant  terms  in 
expressions  for  My  and  Mz  in  the  torsion  equation  but  retain  the  more 
accurate  expression,  including  double-underlined  terms,  for  My  and  Mz  in 
the  bending  equations.  In  the  present  mixed  formulation,  however,  dropping 
of  the  double-underlined  terms  (k2  =  0)  in  Equation  (7)  would  result  in 
matching  the  torsion  equation  of  a  displacement  formulation  equation  but 
result  in  an  oversimplified  set  of  bending  equations.  Retention  of  the 
double-underlined  term  (k2  *  1)  results  in  matching  the  bending  equations 
of  a  displacement  formulation,  but  adds  additional  0(e3)  terms  that  are 
inconsistent  in  the  torsion  equation  of  the  displacement  formulation.  Thus, 
a  consistent  set  of  equations  in  a  mixed  formulation  may,  when  transformed 
to  a  displacement  form,  result  in  neglect  or  retention  of  terms  that  would 
be  differently  treated  in  a  set  of  equations  written  consistently  in 
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TABLE  1.  VALUES  OF  INERTIAL,  ELASTIC,  AND  GEOMETRIC  PARAMETERS  FOR 
COMPARISON  OF  RESULTS  WITH  [1]  AND  [9]. 


(a)  Reference  [1]  Configuration 
p^cR/m  *  5/ (6tt) 
bc/(irR)  “0.1 
c/R  =  tt/40 


a  *  2tt 


Cj  / a  =  0.01 
°o 

EIy/(m07R4)  =  0.014605  (ww  =  1.150) 

10.026787  (wv  =  0.70,  soft  in-plane  case) 
0.16696  (ufy  =  1.50,  stiff  in-plane  case) 
GJ/(m02R4)  =  0.0056732  (u^  =  50) 

VR  =  km2/R  -  0.025 

km1  *  eA  -  e  “  0 

(Vkm)2  -1-5 

0  *  0,  0.3  rad 

(b)  Reference  [9]  Configuration 


0  =  0 
Poo  =  0 

R  =  40  in. 

m  =  0.000125  lb  sec2  in.-2 
EIy  =  25000  lb  in.2 
EIZ  =  7500  lb  in.2 


0  =*  tt/4  rad 

GJ  =  9000  lb  in.2 

km1  “  k,„2  =  1  in. 

eA  "  kA  =  ° 
e  =  </2  in. 


(c)  Tapered  Case 

All  properties  at  blade  root  identical  to  those  listed  above  in  part  (a) 
for  the  soft-in-plane  case.  All  properties  are  constant  except  Ely, 
EIZ,  GJ,  and  m,  which  are  multiplied  by  the  factor  (1  -  O.lx). 


displacement  form.  We  expect  differences  such  as  these  to  be  small. 
Otherwise,  we  must  conclude  that  whatever  has  been  taken  as  an  ordering 
scheme  must  itself  be  somehow  inconsistent. 

A  sample  of  results  generated  for  the  simplified  configuration  from 
[9]  is  presented  in  Table  2.  The  rotor  speed  is  zero  in  this  case,  and 
therefore  there  is  no  aerodynamic  loading  and  no  static  equilibrium  defor¬ 
mation.  There  is  nevertheless,  an  offset  between  the  blade  mass  center  and 
shear  center  and  nonzero  pitch  angle  (the  properties  are  given  in  Table  1). 
Thus,  all  degrees  of  freedom  are  coupled.  The  number  of  segments  used  by 
the  example  from  [9]  is  unknown.  However,  it  is  clear  that  the  present 
results  based  on  a  second-order,  one-step  method  (C'  =  0  for  this  case, 
since  the  beam  is  uniform  with  no  tension)  are  tending  towards  those  of 
[9]  as  the  number  of  segments  is  increased.  With  only  24  segments,  the 
first  mode  is  within  0.1%  of  Murthy's  result. 


TABLE  2.  COMPARISON  WITH  RESULTS  FROM 
[9]  FOR  FREE  VIBRATION  FREQUENCIES  OF 
BLADE  WITH  MASS  CENTER  OFFSET. 


w(rad/sec) 

Segments 

Mode  1 

Mode  2 

Mode  3 

16 

30.7552 

53.6968 

179.8088 

24 

30.7962 

53.7691 

182.3980 

32 

30.8107 

53.7947 

183.3519 

40 

30.8174 

53.8065 

183.8023 

Ref,  [9] 

30.8295 

53.8277 

184.6175 

When  the  rotor  speed  is  nonzero,  the  presence  of  steady  aerodynamic 
loads  and,  in  some  cases,  inertial  loads,  causes  significant  static  defor¬ 
mation  which  must  be  taken  into  account  properly  in  order  to  obtain  correct 
stability  eigenvalues  [1].  Results  from  the  present  analysis  (COLSYS)  are 
presented  in  Table  3  along  with  those  of  [1]  for  both  soft-in-plane  (u>v  =  0.7) 
and  stiff-in-plane  (u^  =  1.5)  cases  with  9=0  and  0.3.  The  quantities 
tabulated  are  tip  deflections  of  v,  w,  $,  and  the  accuracy  is  specified  to 
be  four  significant  figures.  The  solution  is  neither  difficult  to  obtain  - 
nor  particularly  time  consuming.  The  agreement  is  quite  good  regardless  of 
the  choice  of  kl  and  k2.  We  conclude  from  this  that  the  static  equilibrium 
solution  is  not  strongly  affected  by  the  presence  of  the  underlined  terms 
in  Equation  (7)  for  the  limited  range  of  values  presented  here.  Such  is 
not  the  case  for  the  stability  eigenvalues,  however. 

Table  4  presents  the  stability  eigenvalues  based  on  a  second-order, 
one-step  method  with  C*  terms  included  and  16  segments  (converged  to 
3  significant  figures).  Again,  the  agreement  is  very  good,  especially  at 
9=0.  The  presence  of  the  single-underlined  terms  (kl  =1)  has  no  effect 
on  the  results  for  9=0  and  very  little  effect  for  0  =  0.3.  The  double- 
underlined  terms  (k2  =1)  raise  the  torsion  frequency  at  9  =  0.3  from  its 
value  at  9=0  instead  of  lowering  it  as  indicated  in  [1].  The  effects 
of  k2  on  lead-lag  damping  are  minor.  The  difference  in  the  torsion  fre¬ 
quency  due  to  the  k2  terms  is  about  2%.  In  a  displacement  formulation, 
these  terms  of  second  degree  in  bending  curvatures  are  neglected  in  the 
torsion  stiffness  (i.e.,  third-degree  terms,  such  as  4>w"  ,  4>v"w",  etc., 
are  neglected  in  the  second-degree  torsion  equation).  This  is  consistent 
in  displacement  formulations  such  as  [1],  [15],  [18],  and  [19].  It  is  not 
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TABLE  3.  VALUES  OF  STATIC  EQUILIBRIUM 
DISPLACEMENTS  AT  BLADE  TIP.  (kl  =  0, 
k2  =  I  for  COLYS  results) 


v/R 

w/R 

$ 

(a)  uv  = 

0.70,  0  = 

0.0 

COLSYS  -0.002321 

0.0 

0.0 

Ref.  [1]  -0.002326 

0.0 

0.0 

(b)  u>v  = 

o.7o,  e  = 

0.3 

COLSYS  -0.03943 

0.09315 

-0.01332 

Ref.  [1]  -0.03940 

0.09314 

-0.01336 

II 

> 

3 

a 

1.50,  6  = 

0.0 

COLYSY  -0.000522 

0.0 

0.0 

Ref.  [1]  -0.000522 

0.0 

0.0 

(d)  wv  = 

1.50,  0  = 

0.3 

COLSYS  -0.03139 

0.09432 

- 

-0.01164 

Ref.  [1]  -0.03129 

0.09412 

- 

-0.01214 

TABLE  4.  COMPARISON  OF  STABILITY  EIGENVALUE  RESULTS  OBTAINED  BY 
USING  VARIOUS  VALUES  OF  kl  and  k2  (16  segments).  Eigenvalues  are 
given  per  unit  0. 


(a)  0  =  1 

D 

kl 

k2 

Lead-lag 

Flapping 

Torsion 

a>v  * 

0.7S2 

0 

0 

-0.0010 

±0.6929  i 

-0.3235 

±1.07881 

-0.3615 

±4.98091 

0 

1 

-0.0010 

±0.6929  i 

-0.3235 

±1.07881 

-0.3615 

±4. 9809  i 

1 

1 

-0.0010 

±0.6929  i 

-0.3235 

1+ 
►— * 

o 

oo 

oo 

-0.3615 

±4.98091 

Ref. 

[1] 

-0.0011 

±0.7014i 

-0.3245 

±1.07511 

-0.3622 

±4.98751 

U>v  ■ 

1.50 

0 

0 

-0.0011 

±1.4938i 

-0.3235 

±1.0767i 

-0.3617 

±4.98221 

0 

1 

-0.0011 

±1 . 4938i 

-0.3235 

±1 .07671 

-0.3617 

±4.98221 

1 

1 

-0.0011 

±1 .49381 

-0.3235 

±1.07671 

-0.3617 

±4.98221 

Ref. 

[1] 

-0.0011 

±1. 50021 

-0.3246 

±1.07411 

-0.3625 

±4.98881 

(b)  9  =  0 

.3 

UJy  “ 

0.70 

0 

0 

-0.0230 

±0.6823i 

-0.3145 

±1.07491 

-0.3552 

±4.95921 

0 

1 

-0.0228 

±0.68191 

-0.3146 

±1.07461 

-0.3564 

±4.97201 

l 

1 

-0.0228 

±0. 68151 

-0.3147 

±1.07431 

-0.3574 

±4.98361 

Ref. 

[1] 

-0.0249 

±0.6931i 

-0.3117 

±1.06761 

-0.3628 

±4.96371 

Ci)y  m 

1.50 

0 

0 

-0.0641 

±1.61431 

-0.2821 

±0.95471 

-0.3520 

±4.95501 

0 

1 

-0.0665 

±1. 56981 

-0.2792 

±0.95451 

-0.3494 

±5.04711 

l 

1 

-0.0667 

±1 .56631 

-0.2789 

±0.95451 

-0.3491 

±5.05241 

Ref. 

[1] 

-0.0632 

±1.58231 

-0.2821 

±0.95071 

-0.3510 

±4.91501 
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obvious,  however,  that  these  third-degree  terms  are  the  only  ones  that  have 
an  observable  effect  on  the  trends.  The  only  way  to  resolve  that  question 
is  to  examine  results  from  a  consistent  analysis  in  which  0(e3)  terms  are 
neglected  with  respect  to  unity,  which  has  not  yet  been  developed. 

4.2.  Convergence  of  one-step  methods 

We  now  address  the  subject  of  convergence  for  the  one-step  methods 
including  the  C'  terms.  In  Table  5,  results  obtained  for  the  soft-in-plane 
case  at  0  =  0.3  are  tabulated  as  a  function  of  the  number  of  segments. 

The  convergence  is  fairly  rapid  at  first  and  tapers  off  as  the  converged 
value  is  approached.  The  smallest  eigenvalues  are  generally  within  1  % 
for  16  segments. 


TABLE  5.  CONVERGENCE  OF  THE  ONE-STEP  METHOD  BASED  ON  EQUATION  (16) 
FOR  VARIOUS  NUMBERS  OF  SEGMENTS  (wv  =  0.7fi,  0  =0.3,  kl  =0,  k2  =  1) 
Eigenvalues  are  given  per  unit  Q. 


Segments 

Lead 

-lag 

Flapping 

Torsion 

4 

-0.0159 

±0 . 7394i 

-0.3099 

±1.20251 

-0.3365 

±4.8960i 

8 

-0.0234 

±0.6902i 

-0.3108 

±1. 09291 

-0.3523 

±4 . 9548i 

12 

-0.0230 

±0.6840i 

-0.3136 

±1.0791i 

-0.3553 

±4. 9676  i 

16 

-0.0228 

±0.6819i 

-0.3146 

±1.07461 

-0.3564 

±4. 97201 

20 

-0.0227 

±0.6810i 

-0.3151 

±1.07251 

-0.3569 

±4 . 9739i 

24 

-0.0227 

±0.6805i 

-0.3153 

±1. 07151 

-0.3572 

±4 . 9750i 

In  the  one-step  method,  the  C'  terms  may  be  neglected  in  some  appli¬ 
cations.  This  is  certain  to  have  an  adverse  effect  on  convergence  for  a 
rotating  beam,  however.  Even  for  a  beam  with  uniform  properties,  the  per¬ 
turbation  equations  will  have  variable  coefficients  due  to  tension  and 
static  equilibrium  terms.  Thus,  it  is  important  to  study  the  convergence 
of  the  stability  eigenvalues  as  a  function  of  the  number  of  segments  for 
methods  which  do  not  use  the  C'  terms.  In  Table  6,  the  lead-lag  and  flap 
eigenvalues  are  presented  for  three  methods  versus  the  number  of  segments. 
Method  (1)  is  the  complete  second-order  method  with  the  C'  terms  included 
(Eq.  (16)).  Method  (2)  is  the  second  order  method  without  the  C'  terms. 
Method  (3)  is  the  first-order  method  (Eq.  (15)).  There  is  little  differ¬ 
ence  in  the  rate  of  convergence  for  Methods  (2)  and  (3) .  The  inclusion  of 
the  C'  terms,  however,  results  in  about  an  order  of  magnitude  reduction 
in  the  number  of  segments  required  for  convergence.  Hundreds  of  segments 
may  be  needed  for  convergence  when  C'  is  neglected  [23].  It  should  be 
noted  that  [6]  uses  only  the  tension  terms  in  evaluation  of  C' .  While 
that  should  result  in  a  substantial  improvement  in  convergence  over  Meth¬ 
ods  (2)  and  (3),  it  would  still  be  inferior  to  Method  (1)  if  the  beam  has 
portions  of  even  modest  nonuniformities  in  stiffness  or  inertial  properties. 

4.3.  Results  for  a  nonuniform  blade 


We  have  included  results  for  a  hypothetical  nonuniform  blade  in  Table  7. 
All  properties  are  the  same  as  those  used  in  the  configuration  for  compari¬ 
son  with  [1],  except  that  Ely,  EIZ,  GJ,  and  m  are  equal  to  their  values 
for  the  uniform  blade  times  the  factor  (1  -  O.lx).  Convergence  for  the 
nonuniform  blade  is  a  little  slower,  and  all  frequencies  are  increased  from 
the  values  obtained  the  uniform  blade.  The  lead-lag  damping  is  only 
slightly  decreased  by  the  presence  of  taper  for  this  case.  This  small 
effect  is  not  surprising,  since  fundamental  bending-torsion  coupling 
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TABLE  6.  EFFECT  OF  C'izi  TERMS  ON  CONVERGENCE  OF  STABILITY  EIGEN¬ 
VALUES  USING  ONE-STEP  METHODS.  (6  =  0.3,  wv  =  0.7ft,  kl  =  0, 
k2  =  1).  Eigenvalues  are  given  per  unit  ft. 


Segments 

Method  (1) 

Method  (2) 

Method  (3) 

Lead-1 ag 

8 

-0.0234 

±0.6902 i 

-0.0201 

±0. 76991 

-0.0181 

±0 . 8398  i 

16 

-0.0228 

±0.68191 

-0.0216 

±0 . 7308  i 

-0.0203 

±0. 7500i 

24 

-0.0227 

±0.6805i 

-0.0220 

±0. 7152  i 

-0.0211 

±0.72441 

32 

-0.0221 

±0.70681 

-0.0214 

±0.  71241 

48 

-0.0223 

±0.6930i 

-0.0218 

±0 . 7009 i 

64 

-0.0224 

±0.69351 

-0.0220 

±0.6953 i 

80 

-0.0224 

±0.69081 

-0.0222 

±0.6920  i 

96 

-0.0225 

±0.6889i 

-0.0222 

±0. 68991 

Flap 

8 

-0.3108 

±1 . 0929  i 

-0.2849 

±1.13671 

-0.2922 

±1  .14311 

16 

-0.3146 

±1. 07461 

-0.3008 

±  1.10601 

-0.3034 

±1.10211 

24 

-0.3153 

±1.0715i 

-0.3059 

±1 .0944  i 

-0.3075 

±1 .09031. 

32 

-0.3085 

± 1.08841 

-0.3095 

±  1 . 0847  i 

48 

-0.3110 

±1.0821 1 

-0.3116 

±1. 07931 

64 

-0.3122 

±1. 07891 

-0.3127 

±1 . 0767  i 

80 

-0.3129 

±1 . 0770  i 

-0.3133 

±1 .0751  i 

96 

-0.3134 

±1.07571 

-0.3137 

±1.07471 

TABLE  7. 

STABILITY  EIGENVALUES  FOR  A  NONUNIFORM 
0  =  0.3).  Eigenvalues  are  given  per 

BLADE  (kl  -  k2  =  1 
unit  ft. 

Segments 

Lead- lag 

Flap 

Torsion 

16 

-0.0230  ±0.70011 

-0.3138  ±1, 08191 

-0.3564  ±5.07981 

24 

-0.0229  ±0.6987i 

-0.314*  ±1.07901 

-0.3567  ±5.08371 

32 

-0.0228  ±0.6983i 

-0.3146  ±1.07811 

-0.3564  ±5.08511 

parameters  [l]  are  not  changed  significantly.  No  significant  computational 
penalty  in  the  calculation  of  the  nonlinear  static  equilibrium  and  the 
stability  eigenvalues  is  incurred  because  of  nonuniformity. 

5.  Concluding  Remarks 

Nonlinear  equations  for  static  equilibrium  deformation  and  linearized 
perturbation  equations  for  small  motions  about  equilibrium,  from  which 
stability  eigenvalues  can  be  obtained,  are  written  in  a  first-order,  spatial 
derivative  form.  These  equations  are  solved  by  COLSYS  [12]  and  one-step 
integration  techniques,  respectively.  Numerical  results  for  uniform  blades 
obtained  from  the  equations  are  compared  with  published  results  [1]  and  [9] 
and  used  to  study  the  convergence  of  the  one-step  methods.  Results  are 
also  presented  for  a  hypothetical  nonuniform  blade.  In  the  course  of  this 
study  several  conclusions  have  emerged. 
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1.  The  equations  of  the  present  mixed  formulation  are  general  enough 
to  treat  nonuniform  pretwisted  rotor  blades  with  chordwise  offsets  between 
elastic  center,  mass  center,  and  shear  center.  Certain  higher  order  cross- 
section  integrals  are  neglected,  and  symmetry  about  the  major  cross-section 
axis  is  assumed  [15].  Neither  the  equations  nor  the  solution  methods,  in 
their  present  form,  apply  to  the  forward-flight  problem. 

2.  COLSYS  is  used  to  solve  the  nonlinear  static  equilibrium  equations 
with  the  present  mixed  formulation.  The  calculation  of  the  static  equilib¬ 
rium  solution  is  neither  particularly  difficult  nor  time  consuming  even 
though  nonlinearities  and  nonuniformities  are  involved.  Thus,  there  is  no 
need  to  limit  the  static  equilibrium  solution  to  some  linear  or  otherwise 
ad-hoc  estimate. 

3.  There  are  slight  differences  in  the  numerical  results  obtained  from 
mixed  and  displacement  formulations  of  rotating  blades  with  geometric  non¬ 
linearity.  The  differences  stem  from  the  equations  that  relate  moments  to 
derivatives  of  bending  rotations.  A  consistent  ordering  scheme  applied  to 

a  mixed  formulation  may  not  produce  exactly  equivalent  equations  when 
applied  to  a  displacement  formulation.  The  main  difference  in  the  results 
is  in  the  magnitude  of  the  torsion  frequency.  The  difference  is  of  the 
order  of  rotations  squared  with  respect  to  unity  (about  2%) .  The  present 
analysis  yields  the  result  that  torsion  frequency  increases  with  increased 
pitch  angle  in  direct  contrast  to  results  in  [1]. 

4.  One-step,  second-order  integration  methods  appear  to  be  a  viable 
means  of  calculating  the  stability  eigenvalues.  In  order  to  obtain  good 
convergence  it  is  necessary  to  include  the  entire  second -order  term,  for 
rotor  blades,  for  which  tension  force,  elastic,  and  inertial  properties  all 
may  vary  along  the  length  of  the  blade. 

5.  The  terms  involving  the  torsion  moment  M*  in  Equation  (7)  were 
found  to  be  negligibly  small  for  the  limited  range  of  parameters  investi¬ 
gated.  The  double-underlined  terms  do  affect  a  basic  trend,  as  described 
in  conclusion  (3).  This  may  indicate  the  need  for  consistently  incorporat¬ 
ing  terms  of  the  next  higher  order  in  elastic  rotations.  Unfortunately, 
the  only  way  to  ascertain  the  correctness  of  this  assertion  is  to  compare 
results  from  such  an  analysis,  which  does  not  yet  exist,  to  those  of  this 
paper  or  ( 1 ] . 
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APPENDIX  A 


Steady  State  Equations 

The  static  equilibrium  equations  are  given  in  the  following  first- 
order  form  with  the  sign  conventions  and  nomenclature  of  [151  As  in  fl5) 

°f  hi8h^r  °rder*  ®l*  8nd  Bj**  warPin«  rigidity  and  shear 

ea^l0nff8e  *nd  Cl*  ar*  ne8lected-  I"  presenting  these 

equations,  care  has  beeft  taken  to  neglect  terms  that  are  0(e2)  compared 

to  unity.  To  avoid  confusion  with  d/dx,  the  primes  have  been  removed  from 
the  subscripts  of  EIy,  EIZ,  and  Vx. 

Lead- lag  equations; 

Vy  "  -"OaIv  +  e  cos(0  +  5)]  -  Ly  (Ai) 

“  -Vy  +  Vxc  +  mefl2x  cos  (8  +  $)  (A2) 

"  l  (A3) 

V  -  (*il  -  aaSjM,  +  (-ax  +  a„*)My  +  (a2  -  28^)^  +  cos(e  +  *)/Z0 

(A4) 

Flap  equations: 

■  -k  (as) 

■  Vz  -  VXS  -  meft2x  sin(8  +  $)  (a6) 

-  3  (A7) 

3'  "  («s?  -  aifOM*  +  (-a,  -  2a1^)My  +  (ax  -  a%$)Mz  +  eAVx  sin(0  +  $)/Z0 

(A8) 

Torsion  equations: 

M’  -  -V  +  VyB  +  mefl2v  sin(0  +  $)  +  mfl2^  -  ^JcosC©  +  *)sin(0  +  *)  -  M 

(A9) 

*'  "  i/G.lM*  +  My(C  +  as9'  sin  0)  +  Mz(6  -  a60'  cos  0)  -  a70'Vx]  (A10) 
where,  for  brevity, 

G0  ■  GJ  +  asVx  -  ag  M  sin  0  +  a$  M  cos  0 

y  z 


Y«  -  Ely 


z0  -  EIZ  -  EA  eA2 

*i  "  ~I(1/Y0)  -  (1/Z0)]sin  0  cos  0 
«2  "  (sin2  0/Y#)  +  (cos2  0/Zo) 


(All) 
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a,  -  (cos2  e/Y0)  +  (sin2  0/Zo) 

«*“*»”  *2 

a.  -  kA2EI_/Z0  . 

5  A  z  9  concluded) 

*«  “  kA2£A  eA^Zo 

a7  -  (EIz//.0)(kA2  -  EJ/EA) 

The  tern  involving  a7  in  Equation  (A10)  has  been  altered  to  agree  with 
the  recent  results  of  [24]. 

In  the  above  equations,  the  tension  in  the  inextensible  blade  is  given 
by  the  uncoupled  expression 


f R 

Vx  ■  fi2  J  out  dx  (A12) 

*x 

The  aerodynamic  loading  is  derived  in  [1].  The  steady  component 
is 

Ly  -  (p«ac/2)  [v12  -  n2x2(cdj)/a)  -  fixv^e  +  $)  ]  (A13) 

L u  -  (pJoac/2)[-Qxvi  +  fl2x2(0  +  4  +  *)  -  O2xv0  +  (02xcB/2)]  (A14) 

-  0  (A15) 

The  term  i  in  Equation  (A15)  is  defined 

*  -  f  CB’  dx  (A16) 

and  is  represented  in  the  analysis  as  the  following  additional  ordinary 
differential  equation 

♦’  -  CB'  (A17) 

with  the  Initial  boundary  value  specified  by  the  values  of  c  and  3'  at 
x(0).  The  induced  velocity  is  given  by 

vA  -  sgn(0  +  ♦#)nR(iro/8)[^l  +  (12/na)  1 0  +  4.J  -  1)  (A18) 


where 

♦8  ■  4  at  x/R  ■  0.75  (A19) 


and 


0  ■  bc/irR 


(A20) 
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APPENDIX  B 


Linearized  Perturbation  Equations 

The  perturbation  equations  given  in  [15]  have  been  converted  to  the 
linear  ordinary  differential  equations  in  first-order  spatial  derivative 
form  with  the  eigenvalue  X.  Although  not  necessary  for  the  method  of 
solution  finally  adopted,  they  are  presented  linear  in  X  as  well. 

Lead-lag  equations: 

Vy  *  m{Xv*  -  eX$*  sin  9  -  ft2[v  -  e$  sin (9  +  $)) 

+  2fl[Xu  -  e(X;  cos  9  +  XB  sin  9)])  -  Ly  (Bl) 

Hi  -  -V  +  Vj  +  C\L  -  mefi2x$  sin(9  +  $)  +  2meGXv  cos  9  (B2) 

z  y  x  x 

v'  -  C  (B3) 

V  -  (axC  -  +  (-ax  +  a4*)My  +  (a2  -  23^)1^  +  ^(a^  -  a2B) 

+  [a„My  -  2axMz  -  (eA/Z0)sin(9  +  *)VX]$  +  eAVx  cos(9  +  *)/Z0  (B4) 

Flap  equations: 

Vz  *  m(Xw*  +  eX$*  cos  Q)  -  (B5) 

»  $z  -  Vx6  -  BVX  -  mefl2x$  cos(9  +  |)  -  2meOXv  sin  9  (B6) 

w’  -  B  (B7) 

B'  "  (ajC  -  a,B)M  -  (aj  +  2ai4>)M  +  (ai  -  ai,$)M  +  M  (ajC  -  aiB) 

*  x  y  z  x 

-  Ua^  +  a„Mz  -  (eA/Z„)Vx  cos(9  +  ♦)]♦+  sin(9  +  *)/Z„ 

(B8) 


Torsion  equations: 

-  -VzZ  -  CV2  +  VyB  +  SVy  -  me [ Xv*  sin(9  +  $)  -  02v  sin(9  +  *) 
-  n2v  ccs(0  +  $)♦  -  Xw*  cos(9  +  $)  +  2JUu  sin  9]  +  mkm2x|* 


+  »fi2o42  -  cos2e  - 


(B9) 


-  (1/GoHftjj  +  [?  +  a,<0'  +  <F ' ) sin  9]fty  +  My£ 

+  le  -  a,<9'  +  5’)cos  e]Mz  +  mzb) 
where  the  constants  a^  and  G0  are  defined  in  appendix  A. 


(B9) 


(BIO) 
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Equations  (Bl)  through  (BIO) 
using  the  relations 

have  been  linearized  in  terms  of  A 

by 

> 

H 

* 

> 

(Bll) 

w*  *  Aw 

(Bl  2) 

$*  ■  xq> 

(Bl  3) 

Additional  differential  equations 
through  (BIO)  are 

needed  to  supplement  Equations  (Bl) 

4-* 

-  C8'  +  £8' 

(B14) 

V*' 

-  Av’ 

(B15) 

w*' 

=  Aw' 

(B16) 

4,* » 

=  A$' 

(Bl  7) 

The  perturbed  tension  and  displacement  in  the  axial  direction  are 

=  -2mfiAv  (B18) 


u'  -  -li  -  88  -  (eA/Z0)(My  sin  0 

-  M  cos  0)  (B19) 

z 

The  aerodynamic  load  contributions  [1]  are 

l*v  -  (p00ac/2){-nxv1$  -  [2Jlx(cd()/a)  +  (0  +  Dv^Av 

+  [2v1  -  nx(6  +  ^)] Aw}  (B20) 

L  *  (p  ac/2){fJ2x2($  +  i>)  -  nzx(v6  +  8v)  +  (S7zxc8/2) 

W  00 

+  [2ftx(0  +  $)  -  v  ]Av  -  OxAw  +  (3/4)cftxA$  -  (c/4) Aw*} 

(B21) 

-  -(Pooac/2)  (  (cz/8)0xA$]  (B22) 

Thus,  Equations  (Bl)  through  (BIO),  (B14)  through  (B19)  form  the  complete 
set  of  16  coupled  ordinary  first-order  differential  equations  required  to 
calculate  the  stability  eigenvalues  of  a  nonuniform,  pretvisted  blade. 
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